Dynamical mean-field approximation for unitary Fermi gas 
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Dynamical mean-field approximation with explicit pairing is utilized to study the properties of 
a two-component Fermi gas at unitarity. The problem is approximated by the lattice Hubbard 
Hamiltonian, and the continuum limit is realized by diluting the lattice. We have found that at 
zero temperature the predictions of this theory for the energy and the pairing gap agree remarkably 
well with the results of full numerical Monte-Carlo simulations. Investigating the evolution of the 
system with temperature we identify the existence of a second order phase transition associated 
with a jump in the heat capacity and the collapse of the pairing gap. 
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Introduction - The properties of a dilute Fermi gas 
with interparticle distance much larger than the effective 
range r e depend on only two parameters. The scattering 
length, a s which is sufficient to characterize the interac- 
tion and interparticle distance or the Fermi momentum 
kp. When the pair interaction is fine tuned to create 
a two-body bound state with zero energy, the scatter- 
ing length diverges and the Fermi momentum remains 
the only length scale. At this point, commonly referred 
to as unitarity, the system acquires universality, as its 
properties become independent of the nature of its con- 
stituents. Dilute neutron gas where a s is an order of 
magnitude larger than r e is a natural example for such 
system. Unitarity conditions can also be achieved with 
cold Fermi atoms fine tuned near a Fesbach resonance 

;u. 

Regardless of the strength of the interaction, the zero 
temperature ground state of a Fermi gas with attractive 
two-body force is a superconductor. As the attraction in- 
creases the nature of the system changes from a BCS su- 
perconductor at weak coupling to a gas of fermionic pairs 
forming a Bose-Einstein condensation (BEC) at strong 
coupling. At unitarity the system is in between these 
two limits and is neither a BCS superconductor nor a 
BEC. 

In this work we apply the dynamic mean field theory 
(DMFTV introduced over a decade ago by Georges and 
Kotliar [2], to study the properties of unitary Fermi gas. 
In this theory, a lattice problem is mapped into a self- 
consistent embedded impurity problem [3|. This map- 
ping becomes exact in the limit of infinite spatial dimen- 
sions d — ► oo due to the localization of the self-energy 
E(fe,w) -> E(w) 4]. For finite dimensions DMFT is no 
longer exact, yet can be regarded as a useful approxima- 
tion in which a purely local self-energy is assumed, hence 
the name dynamical mean field approximation (DMFA) . 

Applying the DMFA to study the properties of con- 
tinuous Fermi gas, we first construct a lattice version of 



the problem [fl and then seek the continuum limit which 
for finite gas densities corresponds to vanishing lattice 
filling. 

The lattice formulation - The many-body Hamiltonian 
describing a dilute low energy Fermi gas is, 

llTL A — ' J 

(7 

where ip a (x) are the fermionic field operators. To con- 
struct a lattice version of this continuum Hamiltonian 
we represent the configuration space as an L 3 cubic lat- 
tice, where L is the number of sites in each spatial di- 
rection. The time direction is kept continuous. Next, we 
replace the position and momentum variables by the grid 
indices x — » n , p — > ^fc , where n, k are integer vec- 
tors. The grid position and momentum are given by an 
and p/a, where a is the lattice spacing. The fermionic 
fields ip a {x) — > (a)~ 3 / 2 ip na are discretized to obey the 
anti-commutation relations {Vw, V4ier} = 5 aa '8 nn i . The 
corresponding lattice theory is the Hubbard Hamiltonian 

a riri' n 

(2) 

where t — 2 m a ' 2 > an< ^ ^ = Vo/a 3 . D is the hopping oper- 
ator, {Di> a ) n = J2i(^ri+e,,a - + ipn-ei,*), where, 
Bj is a unit vector in the direction j. The spectra of the 
free lattice Hamiltonian is given by 
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In the following we shall use natural units setting h = 
m = 1. The strength of the two-body interaction can 
be related to the two-body scattering length a s through 
summation of the ladder diagrams for two-fermions inter- 
acting at zero energy, zero temperature and zero chemical 
potential, /it — > 0~, [1,0] 

' 1 1 ' (4) 
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where 
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At unitarity Vq c — ~2a/C, so U c = — 7.91576t. In con- 
trast to the original continuum Hamiltonian |T]), the lat- 
tice theory ([2]) acquires an effective range due to the finite 
lattice spacing. The ratio between this effective range, 
r e jf w 2a /n 2 [1], and the average interparticle distance 
is approximately \/n, where n is the average lattice 
filling, i.e. the number of particles per site. For finite 
lattice filling 0.1 > n > 0.01 the corresponding average 
interparticle distance is roughly 12 to 24 times the effec- 
tive range. 

DMFA - Integrating out the fermionic degrees of free- 
dom on all lattice sites but one - the impurity site - the 
DMFT maps the many-body Hamiltonian @ into a sin- 
gle site effective action determined self-consistently from 
a bath with which the impurity site hybridizes. The im- 
purity site can be taken to be a single lattice node Q 
or a cluster of nodes Q. Large cluster size improves the 
accuracy of the DMFA. As the computation complexity 
of the DMFA grows substantially with cluster size it is of 
great interest to asses the quality of the simplest approx- 
imation within this framework, namely the single node 
impurity. 

Using the Nambu formalism, for a system with a super- 
conducting long-range order, the impurity site effective 
action takes the form 0] 
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where [3 = 1/T is the inverse temperature, ^ = (c|,C|) 
are the Nambu spinors, and Go( T ) is given by 



~ ( FS(r) -Qo(-r) 



(7) 



In the following we shall use the "hat" notation for the 
Nambu matrices. The corresponding impurity Green's 
function is given by 

Q{r) ee -<T*,(r)*|(0))s a „ = ( ^ 

'(8) 

where Q{t) = -{Tc„{t)<& (0)) Se// and T{t) = 
-(Tc T (r) Ci (0)) Sef/ . 

In the DMFA the interaction effects are taken into ac- 
count through the self-energy matrix 



E(iw„) 



(9) 



S(iu> n ) -E*(mj„) 
deduced from the Dyson equation 

T,{iuj n ) = QQ 1 {iuj n )-g~ x (iu) n ), (10) 



where Lo n = (2n + 1)tt//3 are the Matsubara frequencies. 
Here and in the following, we have assumed that the sym- 
metry of the pairing is such that the off diagonal self 
energy obeys S(iu> n ) — S*(—iu> n ). The connection to 
the physical lattice is made through the self-consistency 
requirement that the impurity Green's function is equal 
to the local lattice Green's function Q(t) = G(t), or 
g(r) =G(t), and F(t) = F(r), where 
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-iu n + jj, - efe - T,*(iu n ) 



iuj n + fi - e k - S(zw n )| 2 -I- S 2 (iu) n ) 
S{iuj n ) 
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\iuj n + fi - e k ~ E(iw„)| 2 + S 2 (iw n ) 



.(11) 



Solving the Impurity Model - Caffarcl and Krauth [101 ] 
proposed to approximate the effective impurity action 
through an Anderson model 

Ti-And = Ho +Hi 

I, a I, a 

+ Y.vMJ-* + c -* a la ) + Un m , (12) 

l, a 

where the interaction of the fermionic field c a with the 
auxiliary bath fermions ai a generate both the normal 
and abnormal components of the "free" impurity Green's 
function Qo^o- This goal is achieved by choosing the pa- 
rameters of the Anderson model ei,Vi, Di to minimize the 
difference between the Go and Q§ nd over a finite range 
of frequencies \ui n \ < u>n- The "free" Anderson's Green's 
function Q^ nd is calculated by numerical inversion of 7i o • 
When the number n s of fermionic fields, et ; and c, is 
smaller than 6 standard diagonalization methods can be 
used to solve TiAnd- For T — the Lanczos method 
makes a calculation with as many as n s — 10 fermionic 
fields feasible 

Extracting the Physics - Solving the DMFA equa- 
tions yields the local approximation for the self-energy 
Y.(k,iLu n ) S(icj n ), and correspondingly 



G(fe, iuj r , 
and 
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F(k,iu n ) = - 



iuj n + fi - e k - £(iw„)| 2 + S 2 (iuj n ) 

S{iuj n ) 



(13) 



\iu) n + fi — efe - £(iw n )| 2 + S 2 (iuj n ) 



(14) 

Once G and F are known, thermodynamic quantities 
such as the number of particles per site, the energy, and 
the superconducting gap can be calculated through the 
Matsubara sums, 
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E=2qJ2 H e l0+ ( l c J „ + efc + M )G(fe,^„) . (17) 

ak n— — oo 

The number of particles n and the gap Ao can also 
be calculated directly from the Anderson's Hamiltonian 
(|T2|) . If Q = G these results would coincide. However, 
at best Q w G, so this is not always the case. Since 
the Matsubara sum contains explicit dependence on the 
lattice density of state it provides a more reliable esti- 
mate for n and Ao. We shall use the relative difference 
S n = (% ^ nAnd)/ n M between the Matsubara density 
tim and the Anderson's model density n^nd as a mea- 
sure for the quality of the impurity solution. For very 
accurate solution of the impurity model we expect that 
5 n <C 1. When <5„ deviates substantially from zero it is 
an indication that the number of auxiliary fields we have 
used in the solution of © is not sufficient. 

Direct evaluation of the Matsubara sums is imprac- 
tical. Having calculated H(iu/ n ), S(iu) n ) over a finite 
range of frequencies we evaluate the thermodynamic 
observables in the following manner. We set Soo = 
iZe(£(iw/v)), and Soo = S(iu)j^), assuming that at u>n, 
the largest frequency we explicitly consider, the self- 
energy has already acquired its asymptotic value. Us- 
ing these quantities we construct the Green's function 
components Goo , Foo by the appropriate substitutions of 
Soo, Soo in Eqs. 1(15]), (TH)) , The Matsubara sums with 
Goo, Foo are evaluated analytically to obtain n^, Eoo etc. 
Then for the frequencies in the range \oj„ \ < ujn we cal- 
culate the difference between G and Goo to obtain 

1 N 

n = n oo + X! ( G ( k ' iuJ n) ~ Goo(fe,icJ„)) , (18) 

" ak n=—N 

and equivalently for any other thermodynamic observable 
of interest. 

Results - The attractive Hubbard model and the BCS- 
BEC cross were originally studied within the DMFT by 
Keller et. al [ll[ who have established the phase dia- 
gram of the system at quarter filling n = 1/2 using a 
half-ellipse density of states. Later on this work was fol- 
lowed by [H,[ll,[l3| who have studied different aspects of 
the transition in the limit d — ► oo. In Q we have stud- 
ied the continuum limit of the metastable Fermi liquid 
phase in d — 3. Here we extend this study to the ground 
state superconducting phase. The properties of a Fermi 
gas with attractive interaction where studied by differ- 
ent groups using different techniques [HI, [lj| [13, [H, [HI • 
At unitarity, it is customary to present the T = ground 
state energy per particle in the form E/N — £Epc where 
Efg = 0-QEp and Ep is the Fermi energy. The value 
i = 0.44 ± 0.01 was calculated by pj, UM, Gjj], using 
the fixed-node diffusion Monte-Carlo (DMC) method. In 
Fig. Q] we present the energy per particle T = DMFA 
results as a function of the interaction strength in com- 
parison with the DMC results. It can be seen that the 



DMFA results are in agreement with the DMC calcula- 
tions. 
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FIG. 1: (Color online) The T = energy per particle as a 
function of the dimensionless parameter l/(a 3 kF). DMFA at 
lattice filling of n = 0.1 with n 3 =4,5 are denoted by filled 
and empty squares. The circle are the DMC of [l7| . 

Since the continuum corresponds to the limit of van- 
ishing lattice filling we explore in Fig. [2] the dependence 
at unitarity of E/N and the gap function Ao on n. As 
the lattice filling approach zero the accuracy of our cal- 
culation deteriorates since with n s < 6 the bath fermions 
are unable to reproduce the fine details of Qo needed for 
such calculation. To demonstrate this point we attached 
to each data point an error bar, S n (E/N) to the energy 
points and <5„Ao to the gap points, where S n is the rel- 
ative density error discussed above. From the figure we 
can see that as expected the error decrease with the num- 
ber of auxiliary fields. It can also be seen that the error 
grows substantially at the lower most densities. Extrap- 
olating the n s = 6 calculation to the continuum limit we 
obtain £ w 0.44 and A ~ 0.64Ep. These values are 
in close agreement to the Quantum Monte-Carlo (QMC) 
calculations of [H, E3, [3 • 

Turning now to study the behavior of the system at 
finite temperature we present in Fig. [3] the energy per 
particle, the chemical potential and the superconduct- 
ing gap at unitarity as a function of temperature for 
lattice filling n — 0.1. From the figure one can clearly 
identify the superconducting phase transition associated 
with the vanishing of the superconducting gap and the 
discontinuity in the derivatives of fj, and E. Similar qual- 
itative behavior was already identified by Bulgac et. al 
[l8| . The critical temperature T c w O.WEp is somewhat 
smaller than their result T c = 0.23Ep but agrees with 
the value T c = O.lbEp of Burovski et. al [20]. In the su- 
perconducting region, the gap function behaves approx- 
imately as A (T) = A (0)(1 - (T/T C ) K ) and the energy 
as E{T) = NEp G (£ + ((T/T c ) x ). Comparing the DMFA 
results with these formulas we have found that K rs 7, 
and A ~ 5. These values are substantially higher than 
the QMC values n — 1.5, A = 2.5 found in [18j. An exper- 
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FIG. 2: (Color online) The T = energy per particle and 
superconducting gap as a function of the lattice filling n. The 
DMFA calculation with n s = 4 are shown with triangles, n s = 
5 with circles, and n s = 6 with squares. Also shown are the 
QMC results of [H, E3, El] - 



a 0.7 

p 

< 0.6 
0.5 

15 

ft 0.4 

m 

Z 0.3 - 

w 




imental indication for a large A can be found in [2l| were 
the value A = 3.73 was reported, for a unitary Fermi 
gas in a trap. It is interesting to note that in contrast 
with the observation of Toschi et. al [22j our calculations 
seems to indicate a coincidence between the disappear- 
ance of the gap and the discontinuity of the heat capacity. 
Thus Ao provides no indications for a phase transition as- 
sociated with the closure of a pseudo gap above T c . On 
the other hand, in our calculations the gap function Ao 
doesn't vanish at T = T c but, like the other thermody- 
namic variables, change its behavior and drops to zero 
exponentially. 



Conclusions - In this work we have applied the DMFA 
to study the properties of a unitary Fermi gas. We have 
found that the predictions of this theory agree remark- 
ably well with the results of full QMC simulations. The 
DMFA energy per particle actually coincide with the re- 
sults of the fixed node Monte-Carlo simulations, whereas 
the gap function is somewhat higher than previous es- 
timates. At finite temperature the DMFA results agree 
qualitatively with those of although the thermody- 
namic functions exhibit a stronger temperature depen- 
dence. These results indicate that the self-energy of a di- 
lute Fermi gas has only a weak momentum dependence, 
since by construction the DMFA equations yields the lo- 
cal approximation for the self-energy S(fc, iui n ) ~ E(iw n ). 
The momentum dependence of £(iw n ) can be investi- 
gated by extending the impurity site from a single node 
into a cluster of nodes. The DMFA results can be further 
refined using a QMC approach at finite temperature or 
the Lanczos diagonalization method at T = 0. 
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FIG. 3: (Color online) The temperature dependence of 
E,fj,,Ao at unitarity, for lattice filling n = 0.1. The energy 
per particle E/N is shown with squares, the chemical poten- 
tial n is shown with triangles, and the gap function Ao is 
shown with circles. Calculations with n a — 6 are shown with 
filled shapes, calculations with n s — 5 with empty ones. 
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